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," Abstract 

We have revisited the calculation of the neutron electric dipole moment in the 
presence of the CP-violating operators up to dimension five based on the chiral 
perturbation theory. Especially, we focus on the contribution of strangeness content. 
In the calculation, we extract the nucleon matrix elements of scalar- type quark 
operators from the results of the lattice QCD simulations, while those of the dipole- 
type quark-gluon operators are evaluated by using the method of the QCD sum 
rules. As a result, it is found that although the strangeness quantity in nucleon is 
small, the contribution of the chromoelectric dipole moment of strange quark may 
be still sizable, and thus may offer a sensitive probe for the CP-violating interactions 
in physics beyond the Standard Model. 



1 Introduction 



The neutron electric dipole moment (EDM) is one of the physical quantities that are 
extremely sensitive to the CP violation in the high-energy theories beyond the Standard 
Model (SM). The contribution of the SM electroweak interactions to the neutron EDM, 
d n , has been evaluated as d n ~ 10 ^ 2 ' c • cm [T||2] . which is still below the experimental 
limit \d n \ < 2.9 x 10 -26 e • cm [3] by several orders of magnitude. Therefore, the neutron 
EDM provides a clean, background-free probe of the CP-violating interactions in physics 
beyond the SM, as well as the so-called 9 term in the QCD Lagrangian. In order to look 
for the neutron EDM, various experiments have been proposed and put into practice. 
Their sensitivities are expected to be much improved mainly due to the use of ultra cold 
neutrons. For instance, the nEDM collaboration at the Paul-Scherrer Institute (PSI) jl] 
aims at giving a sensitivity of ~ 5 x 10~ 27 e cm, and finally getting into the regime of 
10~ 28 e cm. 

From a theoretical point of view, on the other hand, one needs to interpret the value of 
or the limits on the neutron EDM provided by experiments in terms of the CP violation 
at parton level. To that end, there have been a lot of previous efforts to derive precise 
relations between the CP- violating interactions and the neutron EDM. At the moment, 
only an approach exploiting the QCD sum rules [5l|6] provides a systematic treatment for 
the derivation, which has been conducted in various literature [7HH]. In this approach, 
it is possible to include the contribution of the QCD 9 term, the quark EDMs, and the 
quark chromoelectric dipole moments (CEDMs) on an equal footing. The QCD lattice 
simulations also offer a promising way of evaluating the neutron EDM induced by such CP- 
violating interactions. Although their results are not precise enough at present [15HT9], 
we expect that the lattice simulations eventually will be able to determine it with high 
accuracy. 

An alternative, in fact a traditional, approach to the calculation of the neutron EDM 
is based on the chiral perturbation theory. This method has the advantage of being able 
to include sea quark effects indirectly, while in the QCD sum rule method such effects 
appear only at the higher orders. The contribution of the QCD 9 term to the neutron 
EDM has been evaluated in Refs. [20TT21] . In Ref. [25], the quark CEDM contributions are 
also discussed. While the Peccei-Quinn (PQ) symmetry is often considered in order 
to suppress the 9 term and solve the strong CP problem, it is important to evaluate the 
neutron EDM induced by the higher dimensional operators, such as the quark CEDMs, 
which are sensitive to TeV-scale physics. In the calculation, the effects of the CP- violating 
operators are incorporated into the CP-odd meson-baryon couplings, which are obtained 
by evaluating the nucleon matrix elements of the CP-even quark scalar operators. The 
authors in Ref. [25J extracted the matrix elements from the baryon mass splittings and 
the nucleon sigma terms also evaluated in terms of the chiral perturbation theory, which 
suggested the sea-quark contribution, especially that of strange quark, might be consid- 
erably large. The recent lattice simulations, on the other hand, show that the nucleon 
matrix element of the strange quark scalar operator is actually quite small contrary to 
one's expectation. Since several independent groups have calculated the values with high 
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accuracy and obtained similar results by using the different methods, the results have be- 
come reliable compared with the previous estimates. Moreover, a latest calculation based 
on the covariant baryon chiral perturbation theory in Ref. [27] actually gives a smaller 
value for the strangeness content of nucleon than those in the previous studies. Indeed, 
it is consistent with the lattice results, while its error is much larger than those with the 
lattice simulations. 

This situation stimulates the reevaluation of the neutron EDM based on the chiral 
perturbation theory with the use of up-to-date results for the nucleon matrix elements. 
After the calculation, we will find that although the strangeness quantity in nucleon 
is small, the contribution of the CEDM of strange quark may be sizable, or, even be 
dominant. This result indicates that the strange quark CEDM still may play an important 
role in probing the CP-violating interactions in physics beyond the SM. 

This paper is organized as follows: we first discuss the CP-odd meson-nucleon cou- 
plings in the presence of the effective CP-violating interactions in Sec. [2j Then, in the 
subsequent section, we derive a formulae for the neutron EDM expressed in terms of the 
CP-odd couplings. In Sec. HJ the nucleon matrix elements which we use in calculating 
the CP-odd meson-nucleon couplings are discussed. Scalar contents of quarks in nucleon 
are extracted from the lattice results, while the dipole-type quark-gluon condensates in 
nucleon are evaluated by using the method of the QCD sum rules. The resultant rela- 
tion between the neutron EDM and the CP-violating parameters is presented in Sec. [5j 
Section [6] is devoted to conclusion. 

2 CP-odd meson-nucleon couplings 

We first write down the effective CP-violating interactions at the scale of 1 GeV which 
consist of the flavor-diagonal operators of light quarks and gluon up to dimension five in 
QCD and induce the CP-odd meson-nucleon couplings: 

C qp = - m q qi9 ql5 q + e G ^-G^G A ^- 1 - d g qg s (G ■ a) l5 q. (1) 

q=u,d,s q=u,d,s 

Here, m q are the quark masses, G A is gluon field strength tensor, G ■ a = G A u a pu T A , 
and G A V = \^^upaG Apa with e 0123 = +1. T A and g s are the generators and the coupling 
constant (a s = g 2 /47r) of the SU(3)c, respectively. The second term of the above ex- 
pression is what is called the effective QCD 9 term, while the third term represents the 
chromoelectric dipole moments (CEDMs) for light quarks. They are dimension-five oper- 
ators, and thus quite sensitive to the TeV-scale physics beyond the SM. The CEDMs of 
light quarks are not only directly generated by the CP-violating interactions in the high- 
energy physics, but also induced radiatively through the integration of the CP-violating 
four-quark operators which include heavy quarks [28]. The coefficients of the CP- violating 
operators, 9 q , 9g, and d q , are all assumed to be quite small, and we keep only the terms 
up to the first order of these parameters. 
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By using the chiral U(1)a transformation, it is always possible to rotate out the QCD 
9 term into the first term in Eq. ([T]). In this work, we exploit this basis, where the effective 
Lagrangian is given as 



= - £ m q qi9' ql5 q 

q=u,d,s 



« ^2 d q qg s (G ■ (7)759. 

q=u,d,s 



(2) 



with ^2 q 6' q = 9 = Oq + ^2 q 9 q . From now on, we only use this basis and omit the prime, 

We still have some degrees of freedom in the choice of 9 q since the SU(3)a chiral 
rotation transforms a set of 9 q into another set. By using the degree of freedom, we take 
a basis such that the tad-pole diagrams of the pseudo-scalar mesons should vanish [20J: 



(n qP |£ qP |$ A ) = 0, ($ a = tt, K, v), 



(3) 



where |Oqp) is the vacuum state in the presence of the CP- violating background sources. 
By using the partially conserved axial-vector current (PCAC) relations, Eq. leads to 
the following conditions: 
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m u m d + m d m s + m u m s 
Also, we parametrize the condensate (qg s (Ga)q) as 



(qg s (Ga)q) = -m (qq) . 



(4) 
(5) 

(6) 



Next, we examine the effects of the CP- violating interactions on the couplings of 
baryons with the pseudo-scalar mesons. The couplings are read off from the pion-baryon 
scattering amplitude caused by the CP- violating operators in the low- momentum limit. 
In Ref. [30], the scattering process accompanied by the creation of an extra pion through 
the CP-violating interaction is also considered. The scattering amplitude of the process, 
however, vanishes in our calculation thanks to the vacuum alignment condition in Eq. ([3]). 
Then, the amplitude is again evaluated by using the PCAC relations as follows: 



S^IAjpIS ') 
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d :i x(B B \[J 5 A0 (0,x),C^(0)]\B C ) , 



C\ 



(7) 
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where f n ~ 92.2 MeV [3TJ is the pion decay constant and J^(x°,x) = Q^^^,T A Q is 
the SU(3)a quark axial vector current with Q = (u, d, s) T the quark triplet. The states 
denoted by Q A and B A represent the meson and baryon octets, respectively. We also 
define the 4x4 matrices A q in the spinor basis as 

1 ~ 

A q = -m q e q t - -d q g s (G ■ a) , (8) 
and write Cqp in the following form: 

£pp = QM-ysQ , (9) 

with 

A = diag(A u ,A d ,A s ), (10) 

being a 3 x 3 matrix in the flavor basis. Then, by substituting it into Eq. ([7]), we obtain 
the scattering amplitude as 

($ A B B \C<p\B c ) = ^(B B \Q{T A ,A}Q\B C ) . (11) 

Now all we have to do is to evaluate the baryon matrix elements in the right-hand side 
of Eq. ( fTTjl . Note that the matrix elements consist of the CP-even operators, though the 
interaction we are interested in here is induced by the CP-violating effects. 

The baryon matrix elements in Eq. (fTTj) are expressed in terms of nucleon matrix 
elements by using the group-theoretical arguments. A detail discussion is given in Ap- 
pendix |A] For convenience, we define the following proton matrix elements: 

S q = {p\qq\p} , D q ee (p\qg s G ■ aq\p) . (12) 

Through this paper, the proton state \p(k)) is normalized as 

(p(k')\p(k)) = (2nf5 3 (k'-k). (13) 

Then, using the equations presented in Appendix [AJ we readily obtain the relations be- 
tween the CP-odd baryon-meson couplings and the CP-violating parameters. Among 
them, we just extract the couplings which include the neutron field n and charged mesons, 
since only such kind of interactions contribute to the neutron EDM. They are given as3 

£j> } = 9pnnPri-K + + g EnK E-nK' + h.c. , (14) 



with 



ipnir 



3y/2f* 



{2X X + X 8 ) 



■ [(m u 9 u + m d 9 d )(S u - S d ) + ]-{d u + d d )(D u - D d 



w = ^^(4Yi + 3y 3 -n) 

' [(m u e u + m s e s )(S s -S d ) + l(d u + d s )(D s -D d )}, (15) 



1 The resultant expressions are consistent with those in Ref. 
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where X t and Y{ are defined in Appendix. |A1 Thus, evaluation of the CP-odd baryon- 
meson couplings reduces to that of the proton matrix elements, i.e., S q and D q . 

Of particular interest is the case where the PQ symmetry [26] is imposed. In such a 
case, the expressions presented above are modified. Since there exist other CP-violating 
sources than the QCD 9 term, 9 is not completely erased by the PQ mechanism but 
effectively induced as [52"|I35] 

^ - t £ k ■ < 16 > 

q=u,d,s H 

Then, the CP-odd baryon-meson couplings lead to 

1 r fflQ ,„ „ s 1 



g\ PQ = -—-[-^(S u - S d ) + -(D u - D d )](d u + d, 



1 r^n ,„ „ , 1 



d 



WIpq = -^ri^iSs - S d ) + -(D s - D d )}(d u + d s ). (17) 

Note that in both Eqs. (115 1 ) and (I17|). the contribution of the strange quark CEDM, d s , 
does not necessarily vanislo even if the strange quark content in nucleon is quite small, 
i.e., S s , D s — > 0. This observation allows us to expect that the contribution of strangeness 
to the neutron EDM is sizable even in such a case, and it will be actually shown in the 
following discussion. 



3 Neutron Electric Dipole Moment 

Before evaluating the proton matrix elements, we deduce a formula for the neutron EDM 
induced by the CP-odd baryon-meson couplings given in the previous section. To that 
end, we first obtain the CP-even vertices which we will exploit to calculate the neutron 
EDM. Such interactions are included in the following effective Lagrangian: 

f 2 

J 7T r 



4 -Tr[D M f/( J D^t/) t ] + Tr[B{ip - M B )B] 
^Tr(57^ 5 {^,5}) - ^TrfBf^S]): 



(18) 



2 v ' ,JL "" J/ 2 

with appropriate mass terms for the meson/baryon mass spectrum. Here, U is defined as 
U = £ 2 with 

1 -7= 



£ = exp 



and 



$ = V / 2$ A T A 



1 Vz^~ Ve 



7T 



7T 



77" 



K+ 
K° 



-2 



(19) 



(20) 



2 On the contrary, the strange CEDM contribution to the isospin-conserving nucleon-?? coupling is 
suppressed in the absence of the strangeness content in nucleon. 
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B = ^^B a X a 



The baryon matrix field B is defined by 

(21) 

with A a the Gell-Mann matrices. The covariant derivatives in Eg. (I18j) are given as 

D^U = d^U - ir^U + iUl„ , 

DpB^dpB + pvB], (22) 



/si _i_ Ai y+ „ \ 



where 

r„ = - + - *ye f L (23) 

is called the chiral connection, and and 1^ denote the external sources for the right- 
and left-handed currents, respectively Further, the definition of £ M is 

^ = i[^(d, - irjt - £(<9 M - UJ?], (24) 

which is referred to as the chiral vielbein. The low-energy constants D and F are 
determined by fitting the semi-leptonic decays B — >■ B' + e~ + z/ e at tree level [54"] : 
D = 0.80(1), F = 0.47(1). From Eq. ( ITS]) , we extract the CP-even meson-baryon inter- 
actions. They are expressed as 

Cbb = ^{d ABC D - tf ABC F)d^ A B B Yl 5 B c } (25) 

with /abc the structure constant of SU (3) and cIabc defined by (Iabc = 2Tr[{T A , T B }T C }. 
Among the terms, 

4"p even = + J^TM^" + _=-(£> - F)^ 75 n^K" + h.c. (26) 

contributes to the neutron EDM. Further, by setting = = —eQA^ with e the electric 
charge of positron, i.e., e > 0, and Q = diag(2/3, —1/3, —1/3), we obtain the interactions 
of photon with mesons and baryons. The relevant terms are as follows: 

Cy = - ieA^n+Ti- - tt+^V" + d^K + R- - K+d^K') 
- e(pAp + - 

+ -^-{D + F)pA l5 nn + - -^-(D - F)E~ A l5 nK~ + h.c. (27) 

V 2/,r V IJtt 

Here, A^ denotes the electromagnetic field. 

By using the interactions obtained above, we now evaluate the neutron EDM. It is 
induced by the diagrams displayed in Fig. [TJ Here, solid, dashed, and wavy lines represent 
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Figure 1: Diagrams which give rise to the neutron EDM. Solid, dashed, and wavy lines 
represent baryons, mesons, and photon, respectively. Further, dots indicate the CP-odd 
interactions given in Eqs. (1141) and (TT5]) . 



baryons, mesons, and photon, respectively. Further, dots indicate the CP-odd interactions 
given in Eqs. ffT4"j) and ffl~5l) . The resultant formulao for the neutron EDM is as follows: 



(D + F)g pnn log ( A ) - (D - F)g SnK log ( 



(28) 



with m n and itik the masses of pion and kaon, respectively. Here, we extract the terms 
which contain logarithmic singularities in the chiral limit. A ~ Mb is a UV cutoff. Fur- 
thermore, there exists tree-level, higher order contribution which acts as counterterms |21j . 
With the contribution added to the above equation, we obtain scale-independent formula 
for the neutron EDM. However, the low-energy constants for the operators which give rise 
to the tree-level contribution are not determined through the symmetry arguments. Thus, 
these unknown factors as well as the lack of other non-singular terms in Eq. (128]) result in 
theoretical uncertainty. In the following calculation, we just directly use the expression 
fT28l) with A taken to be the proton mass as an approximation. The error originating from 
the approximation is to be estimated by varying the cutoff scale, similar to Ref. [2T] . 

4 Evaluation of S q and D q 

Now all we have to do is to evaluate the proton matrix elements, S q and D q . 

4.1 S q 

As we mentioned in the Introduction, S q is directly extracted from the lattice QCD 
simulations. To evaluate them, we first introduce the following parameters: 

m u + m d - 
VnN = {p\uu + dd\p) , 

a s = m s (p\ss\p) , 
^ = (P\uu - dd\p) 
(p\uu + dd\p) 



3 



Note that the formula presented here is different from that in Ref. by a factor of two. 
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Exploiting these parameters, we express S q as 



Su = o7 , 7 (1 + Qv-kN, 

2 (m u + m d ) 



1 2 

3d = o7 ; r(l _ a *N, 

2 (m u + m d ) 



= — • (30) 



Recent lattice simulations predict 

a nN = 45 ± 6 MeV, 
a s = 21 ± 6 MeV. (31) 

The symmetry breaking parameter £ is obtained from the baryon mass splittings as [37] 

_ m H o + m 5 - - m s + - m s - / m u + m d 2a s \ 

m s o + 777,3- + m s+ + ms- — 2m p — 2m n \ 2m s a^N J 

( m,, + m,i 2o„ \ , . 

= 0.197 1 — - . (32) 



2m s cr nN 

Then, with the quark-mass parameters [31 
m u + m d 



3.8 ± 0.6 MeV, m s = 95 ± 5 MeV, (33) 



we obtain the values of S q as 



S n = 7.0, S d = 4.8, S s = 0.22. (34) 



Note that the values of S u and S d are almost as twice as those in Ref . [25] . This is mainly 
due to the quark- mass parameters. Moreover, S s is much smaller than that in Ref. [25] . 
since recent lattice simulations predict relatively small values for a s . This fact drives us 
to reevaluate the strange quark contribution to the neutron EDM, as we mentioned to in 
the Introduction. 



4.2 D q 

For the values of D q , on the other hand, we are not able to use the lattice simulations at 
the present time. Although it is desirable that the lattice simulations eventually determine 
D q , for the time being, we evaluate them by using the method of QCD sum rules [HIE]. 
A similar approach is studied in Ref. [30] . 

For this purpose, we first add to the ordinary QCD Lagrangian the following terms 
which couple to the external c-number fields J q : 

£ ext = J q Q9sG ■ aq . (35) 
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The method of QCD sum rules requires that we evaluate the correlation function of the 
proton interpolating fields in terms of two different ways; one is to describe it from the 
phenomenological point of view, and the other is to calculate it in terms of the operator 
product expansion (OPE). The correlation function we are to evaluate is defined as follows: 

IL(k)\j = i J d 4 xe ik - x (n\T{ Vp (x),fJ p (0)}\n)j , (36) 

with rj p (x) the interpolating field for proton. The subscript of the correlator, J, indicates 
that the function in Eq. ( 13 6 p is evaluated in the presence of external fields. We expand 
the correlation function in terms of J q as 

u{k)\j = n°(k) + j q ii q {k) + o(j 2 ), (37) 

and focus on the correlator U q since it includes information on D q . 

The matrix element of the proton interpolating field between the vacuum and the 
one-particle proton state is given as 

(n\r, p (x)\p(k, s)) = ~^=u p {k, s)e- ik - x . (38) 

Here, u p (k, s) is a spinor wave function, which is normalized as usual, i.e., u(k, s')u(k, s) = 
2m p 5 ss '. We parametrize the field renormalization constant as X p , whose value is to be 
determined later. 

A phenomenological description of the correlator is readily given as 

U q{ k) = XlD, f ±- J ±^- + ..., (39) 

where the dots indicate the contribution of the excited states. Among the terms in the 
above expression, we hereafter focus on the terms that contain only one gamma matrix, 
that is, the terms proportional to ^. Then, it follows that 



Il q (k) = ^f q {k 2 ) + ... , (40) 



with 



^^tM^,,,,, (41) 

[k z — m p Y k z — m p 

Here, A q (k 2 ) and B q (k 2 ) are the functions that have no pole at k 2 = m 2 . Dots in Eq. ([4*0]) 
correspond to terms with no gamma matrices. 

In order to evaluate the correlation function by using the OPE, on the other hand, 
one needs to express the proton interpolating field in terms of a composite operator of 
quarks which has the same quantum numbers as those of a proton. The general form of 
the proton interpolating field is given as 

Vp(x) = jpi(x) + /3j p2 (x) , (42) 
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where 

j pl {x) = 2e abc [u T a {x)C lb d b {x)}u c {x) ) 

j p2 (x) = 2e abc [u T a (x)Cd b (x)} l5 u c (x). (43) 

Here the subscripts, a, b, c, are the color indices and C denotes the charge conjugation 
matrix. Although the interpolator j P 2{x) vanishes in the non-relativistic limit, there is 
no reason for the interpolator to be excluded since light quarks in a proton are actually 
relativistic. The unphysical parameter /3 is to be fixed later. The OPE calculation is 
conducted with the quark propagators and condensates in the presence of the external 
fields. They are derived in Appendix [Bj By using them, we evaluate the correlation 
function and extract the terms proportional to the external fields. First, we decompose 
the correlation function as 

n(x)| J = <n|r{» fr ( a :)^(o)}|n> J 

= (Q\T{ Jpl (x)] pl (0)}\Q)j + 0{(n\T{j p2 (x)j pl (O)}\n)j 

+ (Q\T{ Jpl (x)] p2 (0)}\Q)j} + p 2 (Q\T{ Jp2 (x)] p2 (0)}\Q)j . (44) 

For convenience, we use the following abbreviation: 

U kl (x) = (tt\T{ 3pk (x)j pl (Q)}\n) , (k, I — 1, 2) . (45) 

They are expressed in terms of the propagators S^ b (x) given in Eq. fl67j) in Appendix. [B] 
as follows: 

n u (x) = 4e abc e a , b , c , [Tr{T m ,(x)j 5 S^)l 5 }S:A^ + S^^Ux^S^)] , 
U 12 (x) = Ae abc e alUd [Tr{S u aa ,(x)^S^ bl (x)}S^ c ,(x)-f 5 + ^ a ,(x)^,(x)7 5 S c ",(x)7 5 ] , 

n 21 (x) = ^^[^{t^ 

n 22 (x) = Ae abc e alb , c , [Tr{^ a ,(a;)^(x)}75^(x)7 5 + l5 S u aal ^)S d bb ^)S^H] , (46) 
where 

S q ab (x) ee CSg(x)&. (47) 
Then, a series of Wick contraction leads to 





» 


\j= 






ni 2 i 




\j = 






n 2 ii 




\j = 




( X 2) 3 ( 3J « + J d)ml 


n 22 i 




\j = 


sJ qq) - 


( 2 n 3 (9J"« + 5J d )mQ, 
(x ) 
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up to the leading order. Here, we keep only the terms including an external field. It is 
followed that 

n (*)U = -^-M)j^Y m l [W 2 + 6 f 3 + 9 ) J « + ( 5 ^ 2 + 2 £ + 5 )M > ( 49 ) 

and its Fourier transformation results in 

n(k)\j = -^i(^)^ lo g(^) + W + 9 K + + 2/5 + 5 ) J J > ( 5 °) 

with A a certain ultraviolet mass scale. 

Then, by applying the Borel transformation to Eqs. ( I40p and (150j) . we now obtain the 
sum rules for D q as 

and _D S = up to the leading order calculation. When we drive the sum rules, we assume 
that the single-pole contributions A q scarcely depend on k 2 around k 2 = m 2 , and regard 
them as constants. Further, we just neglect B q with expecting that their contributions 
are sufficiently reduced by the Borel transformation. Note that the sum rules in Eq. (l5Tj) 
contain the identical function of M 2 , M 4 exp(mp/M 2 ), in their right-hand side. Moreover, 
the tangent line to the function at a given Borel mass squared M 2 gives us both the first 
and second terms in the left-hand side of the sum rules [14J. Especially, if one sets the 
Borel mass at the minimal of the function, i.e., M 2 = m 2 /2, the single-pole contributions 
vanish. Since this choice also makes the sum rules stable under the slight variation of the 
Borel mass, we adopt it in the following discussion. 

Now all we have to do reduces to the determination of the constant A p . One way to 
evaluate it is again using the method of QCD sum rules. In Ref. [38], two sum rules which 
allow us to extract A p are presented; one contains the Lorentz structure j$ and the other 
is proportional to unity. On the other hand, it is possible to extract the value from the 
lattice result. By using the parameters obtained in Ref. we obtain 

A p = 2 x 0.91 • [(«! - a 2 ) - p( ai + a 2 )] , (52) 

with the parameters a\ and a 2 given as 

ai = -0.0112 ± 0.0012 (stat) ± 0.0022 (syst) GeV 3 , 

« 2 = 0.0120 ± 0.0013 (stat) ± 0.0023 (syst) GeV 3 . (53) 

Here we take the renormalization effect into account [2]. In the following calculation, we 
use the values for \ p evaluated in both approaches, and present two results corresponding 
to these two values. 
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The estimation of D q presented here is, however, to be considered as just for reference. 
At the moment we only present center values for those quantity. In order to estimate the 
theoretical error for the calculation and to improve the accuracy of the result, we need the 
investigation of the excited/continuum state contribution as well as the execution of the 
higher order calculation for the OPE. In the higher order calculation, there exist several 
unknown condensates such as the susceptibilities of (qq) with respect to the external 
sources, and thus they make it difficult to give a precise prediction. Therefore, much 
detail analysis with the determination of these unknown condensates is required for a 
robust calculation of D q . We hope our study presented in this paper stimulates a lot of 
efforts to evaluate the quantity with various methods, especially with the lattice QCD 
simulations. 

5 Results 

Now we derive a relation between the neutron EDM with the CP-violating parameters by 
using the nucleon matrix elements obtained above. We use the values in Eq. (154"|) for S q , 
while D q (q = u, d) are evaluated from the sum rules in Eq. (15"Tj) with D s presumed to be 
zero. 

We are interested in the significance of the strange CEDM contribution to the neutron 
EDM. In order to examine it, we consider each contribution of CEDM in the presence of 
PQ symmetry. In this case, the resultant relation between the CEDMs and the neutron 
EDM is expressed as 

d n \pQ = y^eC q -^. (54) 
q H 

Here, we normalize the CEDMs by the quark masses, since in most cases the quark 
CEDMs are proportional to the quark masses, d q oc m q . Then, we look into the ratio of 
the coefficient C s against those of up and down quarks. It allows us to evade an overall 
factor of uncertainty coming from the logarithmic factor in Eq. (1251) . though the relative 
values of C q still might be affected. 

In Fig. |2j we plot the ratios C s /C q as functions of the parameter ft in the sum rules. 
The upper blue and the lower orange lines show C s /C u and C s /Cd, respectively. The 
solid lines correspond to the results in which we use X p in Eq. (132]) . while the dashed 
lines represent those obtained with X p evaluated in Ref. [55] . In the calculation, we take 
ml = 0.8 GeV 2 [29] and (qq) is evaluated using the relation (qq) = —m\f1j{m u + m d ) as 
(qq) ~ —(280 MeV) 3 . In addition, when we compute X p from the results in Ref. [55], we 
use the sum rule for j> and set the threshold w and the gluon condensate b = (2ir) 2 (^-GG) 
as w = 2 GeV and b = 470 MeV 4 , respectively. It is found that both calculations for 
X p give rise to similar results in the case of ft < 0, while one deviates form the other 
for positive ft. In any case, the strange CEDM contribution is sizable, in fact dominates 
the other contributions, though the strangeness content of nucleon is quite small. The 
contribution comes from the i^-meson loop and does not vanish in the limit of small S s 
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Figure 2: Ratios C s /C q (q = u, d) as functions of the parameter (3 in the sum rules. Upper 
blue (lower orange) lines show C s /C u (C s /Cd). Solid lines correspond to the results where 
we use \ p in Eq. (I5"2"j) . while dashed lines represent those obtained with \ p evaluated in 
Ref. §E\. 



and D s , as it is understood from Eq. (fT7I) FI 

As we have mentioned to above, the computation suffers from the uncertainty due 
to the tree-level contribution of the unknown counterterms as well as the lack of other 
non-singular terms in Eq. (I28p . To estimate its significance, we vary the cutoff scale A in 
Eq. ( 12"5|) . following the analysis in Ref. [T2]. In Fig. El we plot the dependence of the ratios 
C s /C q on A. The blue solid (orange dashed) line shows the result for C s /C u (C s /Cd). 
We use the lattice result in Eq. fl52|) for the value of X p with (3 set to be —1. This figure 
tells us that although the variation of the cutoff scale may change the ratio by an 0(1) 
factor, the contribution of strange CEDM tends to remain important. Taking the results 
obtained above into consideration, we conclude that although the recent lattice results 
indicate that the amount of strange quark in nucleon is small, the CEDM of strange quark 
may play an important role in probing the CP-violating interactions in physics beyond 
the SM with the neutron EDM. 

Finally, as a reference value, we present the numerical results for the relation between 
the CEDMs and the neutron EDM: 

d n = e(5.8 x 1(T 16 # [cm] - 1.5d u + 0.6d d + 0.44), (55) 

where we take /3 = — 1 and use A p in Eq. (I52p . Note that unlike the CEDM contribution, 
the calculation of the 9 contribution is not affected by the uncertainty in the values of D q , 

4 The strange CEDM contribution to the neutron EDM vanishes if m^Sd + Dd = 0, though it looks 
accidental. 
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Figure 3: Dependence of C s /C q on cutoff scale A. Blue solid and orange dashed lines show 
the results for up and down quarks, respectively. We use \ p with [3 = — 1 in Eq. (|52|) in 
the calculation. 



since the contribution only depends on S q . In this sense, the result for the 9 contribution to 
the neutron EDM is robust up to the uncertainty coming from the chiral loop calculation. 
In the presence of the PQ symmetry, on the other hand, the relation is modified to 

gUpq = e(3.6d u + 3.1d d + 0.5d s ). (56) 

Compared to the results obtained with the QCD sum rules [T2"|[T4"] . our results here predict 
relatively large neutron EDM. Indeed, a latest calculation conducted in Ref. [14] give^l 

d n = e(8.2 x 1(T 17 £ [cm] - 0.37d u + 0.37d d - 0.02d s ), (57) 

and in the presence of the PQ symmetry, it reduces to 

d n | PQ = e(0.36d u + 0.714). (58) 

These values are found to be smaller than ours by 0(1) factors. With the present knowl- 
edge, we are not able to conclude which method is appropriate for the calculation of the 
neutron EDM, since both of them may have substantial uncertainty. 



6 Conclusion 

We have calculated the neutron electric dipole moment in the presence of the CP-violating 
operators up to dimension five based on the chiral perturbation theory. In the calculation, 

5 Equations ([57)1 and ([55)1 also contain the contribution of quark EDMs; 0.95^ — 0.24^ with d u and d<j 
the EDMs of up and down quarks, respectively. Further, the numerical values presented here are in fact 
different from those in Ref. [TJ] by nearly a factor of two. The difference stems from the use of different 
values for the quark condensate; (qq) = —(225 MeV) 3 in Ref. [Tl] . 
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we extract the nucleon matrix elements of scalar- type quark operators, S q , from the recent 
lattice results, while those of the dipole-type quark-gluon operators, D q , are evaluated by 
using the method of the QCD sum rules. Especially, we focus on the strange CEDM 
contribution to the neutron EDM. Although the lattice QCD simulations tell us that the 
strangeness content in nucleon is small, it is found the contribution may not be suppressed. 

There are two origins of the theoretical error in our calculation; one is unknown coun- 
terterms in the effective Lagrangian which contribute to the neutron EDM at tree level, 
and the other is the evaluation of D q by using the method of QCD sum rules. The former 
is somewhat inevitable until one fixes all of the low-energy constants for the counterterms. 
The latter is, on the other hand, expected to be improved if the higher order calculation is 
conducted with unknown condensates determined from other studies. Further, the lattice 
QCD simulations might compute D q directly. It allows us to make a prediction of neutron 
EDM much precisely with error estimated more systematically. 

At present, we conclude that the contribution of strange quark CEDM may be still 
significant, and therefore it might offer a sensitive probe for the CP- violating interactions 
in physics beyond the Standard Model. 
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Appendix 

A Baryon matrix elements 

In this section, we derive formulae for the baryon matrix elements given in Eq. (ITT]) . To 
that end, we first rewrite A in terms of the SU(3) generators: 

A = -{An + A d + A s )l + {A u - A d )T 3 + -±=(A U + A d - 2A S )T 8 . (59) 

Then, by using the anticommuting relations 

{T A ,T B }= l -5 AB + d ABC T c , (60) 
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we obtain 

(B B \Q{T A ,A}Q\B C ) 

= 2 -{B B \Q{A u + A d + A S )T A Q\B C ) + ^8 A3 (B B \Q{A U - A d )Q\B c ) 

+ d A3D (B B \Q(A u - A d )T D Q\B c ) + ^=8 AB {B B \Q(A V + A d - 2A S )Q\B C ) 

+ l=d A8D (B B \Q(A u + A d - 2A S )T D Q\B C ) . (61) 

Here, we define d ABC as d ABC = 2Ti[{T A ,T B }T c ]. 

Next, we express the baryon matrix elements in the above equation in terms of the 
baryon matrix field B defined in Eq. (j2lj) . Since there are two different ways of combin- 
ing three matrices to form an SU(3) covariant form, we parametrize the baryon matrix 
elements as 

(B B \Q(A U + A d + A S )T A Q\B C )WB C = X{Ti(BT A B) + Y{Vr(BBT A ) , 
(B B \Q(Au - A d )T A Q\B c )WB c = X 3 Tr (BT A B) + Y 3 Tr(BBT A ) , 
{B B \Q(A U + A d - 2A S )T A Q\B C )WB C = X s Ti(BT A B) + Y 8 Ti(BBT A ) . (62) 

Further, considering the combination of two matrices to form an SU(3) singlet, we obtain 

(B B \Q(Au - A d )Q\B c )WB c = Z 3 Tv(BB) , 
(B B \Q(Au + A d - 2A s )Q\B c )Wb c = Z 8 Ti(BB) . (63) 

These parameters can be written in terms of proton matrix elements. Then, using S q and 
D q defined in Eq.f lT2"|) . we express the parameters in Eqs. f lo^]) and flo^]) as 



x x = 


-{m u 9 u 


+ m d 9 d + m s 9 s )(S u - Sd) -^(d u + d d + d s )(D u - D d ) , 


x 3 = 




- m d 9 d )(S u - S d ) -~(d u - d d )(D u - D d ) , 


x 8 = 


-{m u 9 u 


+ rn d 9 d - 2m s 9 s )(S u - S d ) - ^(d u + d d - 2d s )(D u - D d 


Y x = 


-(m u 9 u 


+ m d 9 d + m s 9 s )(S s - S d ) - -(d u + d d + d s ){D s - D d ) , 


Y 3 = 


-{m u 9 u 


- m d 9 d )(S s - S d ) - - d d )(D s - D d ) , 


Y s = 


-(m u 9 u 


+ m d 9 d - 2m s 9 s ){S s - S d ) - -{d u + d d - 2d s )(D s - D d 


z 3 = 


-{m u 9 u 


- m d 9 d )(S u + S d + S s ) - ^(d u - d d ){D u + D d + D s ) , 




-{m u 9 u 


+ m d 9 d - 2m s 9 s )(S u + S d + S s ) - ^(d u + d d - 2d s )(D u 



(64) 
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B Quark propagators and correlators of the back- 
ground fields 

In the OPE calculation carried out in Sec. I4.2[ we need to obtain the quark propagators 
as well as the correlators of the quark/gluon background fields in the presence of the 
interaction in Eq. ( 13 5p . The quark propagators are defined as follows: 

[S q ab (x)U = (V\T [q aa {x)q bp m , (65) 

where a and (3 denote spinor indices. Furthermore, we perturbatively expand the propa- 
gators as 

[Sli^U = [ S T(x)] af} + XIMXUO) + [S q ab (x)] a , \ g +... ■ (66) 

The first term is the free propagator, and the second term describes the correlator of 
the quark background fields, with x q aa ( x ) a classical Grassmann field which indicates the 
quark background field. The third term represents the propagator which includes one 
gluon. Let us evaluate these terms in x-space. The first term, [S q } b Q \x)] a p , is readily 
evaluated as 

s m, x) _ ^abji_ (67) 

^ab \ X ) - (2-2)2' 

where we neglect the quark masses since their contribution only appears in the higher 
order operators. 

Next, we evaluate the third term in Eq. ( 1661) . In this calculation, it is convenient to 
exploit the Fock-Schwinger gauge [2] for the gluon field. In this gauge, the gluon field is 
subjected to the following gauge fixing condition: 

x^G^x) = , (68) 

where G^(x) = G A (x)T A is the gluon field. Then, it is expanded by its field strength 
tensor such that 

G^x) = ^/G ¥ (fl) + ^//(I) a ^(0)) + ^AV(l) a ^ )1 (0)) + .... (69) 
By using the expression, the gauge covariant form of the propagators is obtained as follows: 



- 2 {t,G ab -a} + U q ^-?t 
x z \x l Y 



(70) 



with G^ h = G A ^ u T A b . Here we keep only the first order terms with respect to the external 
fields, J q . 

Also, we translate the quark and gluon background fields into their condensates. The 
correlation function of quark background fields, X<L( x )Xb,a(0) > i s related with the quark 
condensate as 

XLWXUO) = (n\T[q aa (x)q bp (0m). (71) 
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By using the Fierz identities and carrying out the short- distance expansion of the quark 
field, 

q ( x ) = q (0) + x^D^O) + . . . , (72) 

we obtain 

xL(x)xUo) = - <gg)[l - JXU ■ ( 73 ) 

In general, the vacuum condensate in the presence of external fields, (qq)\j, is different 
from ordinary one, (qq) [12]. However, the susceptibilities of the condensates to the 
external fields only appear in the higher order of the OPE, therefore, we neglect them 
since in our calculation we just deal with the leading order contributions. 

In addition, we need the interaction part of the quark and gluon background fields, 

9sXla(x)xU0)G£ u = (Q\T[g s q aa (x)G^(0)]\Q), (74) 

and it leads to the following equation: 

9sXL(x)xU0)G^ = -Lr^mKqq^U. (75) 
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